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We obtain the exact nonequilibrium work generating function (NEWGF), for a small system 
consisting of a massive Brownian particle connected to internal and external springs. The external 
work is provided to the system for a finite time interval. The Jarzynski equality (JE), obtained in 
this case directly from the NEWGF, is shown to be valid for the present model, in an exact way 
regardless of the rate of external work. 
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I. INTRODUCTION 

Due to the development of precision manipulation techniques at very small length scales, such as atomic force 
microscopes or optical tweezers, it has become possible to study the response of small systems to applied external 
influences, such as the work done by pulling apart the extremities of DNA molecules [l|-|3(- These nonequilibrium 
experiments can yield equilibrium information (free energies) about the system [3|45|. This can be obtained from 
relations such as the Jarzynski equality [(| @] (JE). The explicit form for the JE Ip relating the external work 
fluctuations and the variation AF of the free energy between these two states is 



exp{-/3W} = exp{-pAF}. (1) 

In order to understand the equality above, we observe that an important step in the derivation of Fluctuation 
Theorems, and the JE, is to make sure that the changes in the values of the Hamiltonian of the system are directly 
associated with the external work performed by the environment onto the system In fact, the distinction between 

environment and system is a matter of choice. We can decide which variables are accounted for as belonging to the 
system, the remaining being part of the environment. In fact, conservative external forces might be incorporated into 
the Hamiltonian of the system as potential-energy terms Q, thus making the total Hamiltonian invariant under its 
effect. So, for an isolated system the external work verifies (x and p represent all phase space variables) 

^external = U( X 2,P2, h) - U(xi,pi, h), (2) 

where positions and momenta evolve from state 1 to state 2 according to the dynamics of the system. However, 
external forces arising from rheonomic constraints (l0| might perform work on the system but, since they do not come 
from a simple expression of a conservative potential energy, they cannot be simply incorporated into the Hamiltonian 
as above. On the other hand, their work expression can be derived from the constraint's equations [l(3]. In fact, these 
forces can do work on an otherwise isolated system and change the Hamiltonian energy landscape upon which the 
phase-space point evolve in time. [4lj Thus, one has to be quite careful when defining the external work, as that choice 
can substantially change the final form for any fluctuation relation Indeed, fluctuation relations that are derived 
under distinct choices of the definition of work used will lead to distinct forms of the fluctuation theorems [1, ||| ITTrl^ 
(FT). 

The JE have been demonstrated exactly for systems initially thermalized at temperature T, that can be either 
mechanically closed or in contact with a thermostat (at temperature T also) all the time. These results are obtained 
for systems that are large enough so that they may be placed into a true equilibrium state for which a valid partition 
function exists, and a correspondingly valid free-energy associated with it (2lj . according to the usual derivation of 
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the canonical ensemble distribution. On the other hand, given the Hamiltonian for a small system, one can define its 
canonical partition function and calculate the corresponding free energy. 

Several models p2T[30j have been proposed where the JE is verified for distinct systems. In particular, for Ref. [22|, 
the system consists of a particle pulled through a thermal bath in a manner that becomes time invariant as t — > oo. We 
propose the present model as a realization of a non-equilibrium process for a system that is clearly non-homogeneous 
in time. Irrespective of the duration of the process, the amount of external work done is finite, in average, and 
fluctuates around a well defined value. 

Our model consists of a Hamiltonian that incorporates the kinetic energy term of a Brownian particle, and the 
potential energy terms associated with two springs connected to the particle: one represents a harmonic potential 
centered at x — an the other spring has one end fixed on the particle while the other end is pulled, according to 
a given time protocol, as work is done on the system (mass and two springs) by the external constraint force at the 
pulling end of the spring. This external force is the only force that can change the energy of the system. Our choice 
of definition for the work follows from the discussion above and seems well suited for studying the JE. This model 
can be thought as a prototype of a controllable heat engine that can operate in reversible or irreversible modes. 

We calculated exactly the model dynamics, in the spirit of other models previously used by the authors [HI [Hj], 
and we analyzed the behavior of external work fluctuations that are in contact with thermal bath. The thermal bath 
is represented by a noise term in a Langcvin equation (LE) [33j | . We approach that problem from the point of view of 
a non-equilibrium work generating function (NEWGF) , equivalent to the complete non-equilibrium work probability 
distribution. Such functions have been used in the context of the JE [34l , |35| . They allow us to obtain the work 
probability distribution for equilibrium and nonequilibrium conditions. 

Starting from the LE for an undcrdamped Brownian particle, we obtain exact information on the structure of the 
NEWGF. Solving exactly the LE for a system where the noise stochastic properties are known is akin to solving 
the Kramers-Moyal equation [33[ for the probability distribution. Establishing exact finite-time Langevin dynamics 
results gives us the possibility to calibrate analytical or numerical models via the JE, in other words, it is a first 
principles calculation that corroborates the validity of JE, and it can also serve as a detailed testing ground for 
numerical simulation models. The model consists of a Brownian particle of mass to, under the action of a harmonic 
potential k, and in contact with a thermal reservoir at temperature T and friction coefficient 7. We attach to that 
particle an external spring (fc'), by one extremity, and pull the other extremity at a fixed time rate (defining the work 
protocol). The particle-reservoir coupling is represented by a Langcvin force (noise) rj(t). The external spring has 
the externally moving extremity at the point x spr i ng (t) — L(t), as work is externally done into the system (to, k, k') 
without ambiguity: the work is the product of the externally varying force applied to the moving extremity of the 
spring k' with its displacement dL. The model and protocol we use, varying an external coordinate according to a 
pre-established time-rate, are equivalent to others found in the literature [9(. 

This paper is organized as follows. In Sec. II we define the model. In Sec. Ill we obtain the generating function for 
the work probability. In Sec. IV we derive the JE, followed by our conclusions in Sec.V. 



II. MODEL 

So, let us define our model LE: 



to v = — 7 v — kx — k (x — L) + r), (3) 
x = v, (4) 

L = L (l — e-^ x ) . (5) 



The process starts (t = 0) with the system initially at equilibrium with a reservoir (at temperature T). The initial 
conditions (xo,vq) are distributed obeying the canonical distribution at temperature T, and with x spr i ng (t = 0) = 0. 
Then, the external spring is moved [given x spr i ng (t) = L(t)} up to t = r, which may or may not be extended to infinity. 
The specific form of L(t) given in Eq.([5]) was chosen for being easy to manipulate, but it can be readily generalized. 
The rate A can be set to any value, with A — > 00 corresponding to the reversible work process. 

By taking the Laplace transform of the Gaussian noise function's, we obtain for the second cumulant (given that 
the average of 77 is null) 

(»K*i)tK*2)) - (6) 
z\ + z 2 



We can integrate this system exactly by techniques similar to the ones used previously j31l . l32l . l36l ]. where the integration 
paths arc all described therein. However, at present, we will use a direct solution technique which is different, and 



simpler than that in references [3l|, |32|, |36[ . 
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III. GENERATING FUNCTION FOR THE EXTERNAL WORK 



In the spirit of the JE we define the precursor function to the NEWGF as 

(-iu) 



F(u) = e X p{-iuW T } = Y^^j^W?, (7) 



n=0 



where the average is taken over the thermal noise, which corresponds to all possible nonequilibrium paths for the 
Brownian particle. In order to construct the NEWGF we need to take averages, thermal (represented by (/)) and 
over the initial conditions (represented by /). The expression for it reads 



T(u) = (F(u) 

£ ^^-({AU + AW p + AW h + AW v ) n ) 

n=0 

exp {-iu (AU + A W p )} £ i^- (A W h ) n £ ((A W v 



where the partial terms for the external work A W„, A Wh, and A W„ will be defined in the following. 

The accumulated work- function, that measures the total external work done on the system up to time r. W T , is 

fet = -#(*(*) --&(*))) 



Wt = J F ext dL = -k'J dt — (x(t) - L(t)) 



= AU-k' [ T dt^-x{t), (8) 
Jo dt 

with A U = fc'L 2 (r)/2. It is the coupling of x(t) and d Q l 'p will give rise to the irreversible work loss. 

A few thermodynamic properties for our system can be obtained directly from the equilibrium partition function Z = 

JZo *)r e~^\ where P H{x t p) = (£ + k x* + k' (x - Lf) . We find F = £ - Tin ' 

S = In {^-j^- ^-jrfp-^j +1, and E = T + ^ £+ k , ^ where T corresponds to the kinetic and elastic energy contributions 

around equilibrium (via equipartition theorem 2 x T/2). The second term on the RHS is the rest-energy of two springs, 
k and fc', of zero length, connected serially with total extension L. 

Keeping T constant, the reversible work W r and the free-energy AF are identical (L(t = 0) = 0): 

equilibrium position x eq = , ^ L(t): 



, 2 , 

The dissipativc work W ( i = W T — W r can be expressed as the integral of fluctuations of x around the instantaneous 

k' 
k+k' 



If) J . (10) 

We can obtain a particular solution xs(t) associated with the source term S(t), via the Green's function method, 
explicitly as 

xs{t) = L dt W^=F sm ( 2 J 5(< >• 

where = 7/m, and w 2 = (fc + k')/m. 

The exact solution x(£) for Eqs. taken in the sense that the exact form for the Langevin term r){t) is known 
(due to the linearity of the problem under scrutiny) is simply the sum of the homogeneous solution Xh{t), the term 
x p (t) with source Lq (l — e -t / A ), and the term x^i) with source rj{t): x(t) = Xh(t) + x p (t) + x v (t). 
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We may write the total work as 

W T = AU-k'f dt^^-(x h (t) + x p (t) + x n (t)) 

= AU + A W£ + A Wp + A . (12) 

The work expressions, for A Wfi° and A W£° are given below. The detailed form of the coefficients is given in 
Appendix E] in Eqs. ED and [MJ 

AWl = Clx + C^v , (13) 
AW; = -k> £ dt?j&x v (t). (14) 
Making the summation for the NEWGF: 



n. 

n=0 k ) 

where xi = and v 2 . = — . 

The thermal contribution for the work arises from the need to compensate the dissipative coupling of the noise 
term with the pulling rate expressed below: 

AW; = -k' j\t^-x v {t). 

The contribution above is the only source of irreversibility into the system. The expression for A is then: 



AW T V = -lim r ^lfj(i qi + e) x 



2L k'\ f T , f* J , e^^-i- 1 ^ . (y/±u*-P(t-t) 
at I at , sm 



^ m J J a Jo \J\tJ 1 - Q 2 

/OO 1 
£ fj(iq + e) (W^(iq + e) + W^(iq + e) + W^q + e) + W^(iq + ej) , (15) 



where the detailed expressions for the WL listed in Appendix lAl 

The Gaussian property of the noise function leads to the disappearing of the odd momenta of AWJ.. The even 
ones are given as products of (A A W 7 ^), that can be expressed as sums of terms of the form: 

hj = lim J™ ^ J™ (f)(i qi + e) ^q 2 + £)) W^iQi + <0 W^(iq 2 + e) 



f°° dq 2 2 7 T 
J-oo 2tt J_ 00 2tt (i qi + iq 2 + 2e) 



% / ^ / ^ - , L , W^(m + <0 W£(«ga + e). (16) 



There are ten possible distinct pairs above but we can show that only In, J 12 , and In give non-zero results. Thus, 
wc have 

{{AW;f n ) = M(/ 11 + 2J 12 + / aa )", 

yielding 



n=0 

The expression for the NEWGF reads 



2^^^(( A ^) > = expj- u \ 



2 C[ 2 ^ + CJ 2 i; 2 2 7 11 +27 12 +7 : 



F(u) = cxp I -iu (AU + A W p )~u 2 1 " 2 2 -u 2 11 ' ^ ' Zl = cxp - M 2 ^ 2 } . 



Generating functions of Gaussian shape [37H39I ]. for quadratic time- varying potentials, have already been found in the 
literature. The explicit expressions for 7?-i and 7^.2 can be found in Appendix [AJ 



5 



IV. THE JARZYNSKI EQUALITY 



The Jarzynski equality (JE) is verified, after some tedious, but straightforward, manipulations of the terms IZi and 
7?. 2 from Appendix for all values of r and A since for u = —i/T, 

<4)-H4))-~(-^fO-->(-")- 

where AF = ■jj^jp^L'j, and Lf = L(t). Despite the highly non trivial dependence of 1Z%, and IZ2, on A and r, at 

u = —i/T the correct cancellations occur and the JE is verified. 

We notice that by fixing the final position of the external spring < Lf < Lq the ratio r/A gets fixed, 

T =h/ Lo 



A \Lo-L 

which gives us an infinite number of distinct protocols for taking the system from state A = (L = 0) to state B = 
(L = Lf), verifying the JE for all cases for a fixed ^f- = Fb 7 f Fa . 

The Gaussian form of F(u) is the expected one due to the linear form of the harmonic potential [4(| • The present 
model is an explicit dynamic solution that could be extended to other forms of the noise, in the case that its cumulants 
are known. In fact, for the quasi-static case, A, r — > 00 (r/A fixed), T{u) — > cxp{— iu AF} — > p(W) = S(W — AF). 
In this case, there is only one way of carrying out the process. 

We can also obtain the forward ratio for the work distributions. We notice that W is finite (with probability equals 
to 1) and positive. The probability distribution for the total work p(W) is the inverse Fourier transform of F\y(u): 

P(W)= I ^F(u)e™ w . (18) 



The expression for p(W) explicitly reads: 



P (W) = W^cxp(-&^lU. 



K 2 I 4^ 2 

It can be seen that the average work done externally on the system is given TZ\, while the variance is given by 27?-2, 
which is proportional to T for all A and r. 

The ratio p(W) /p(—W) , not to be mistaken with the Crooks Fluctuation expression Q, can be obtained explicitly 

p(W) fWUi 

"2 



P (- w) -«n«r J - (20) 



The expression above is well behaved in the instantaneous work case, A — > 0, since averaging over an ensemble of 
initial conditions distributed with temperature T converges, while it becomes proportional to a delta-function when 
A — > 00, since 7?-2 — > in this case. In fact, the averages of the work do depend on the work rate A and the elapsed 
time t, but the JE arises regardless of it. 



V. CONCLUSIONS 



In conclusion, we develop an exact technique appropriate for treating a system consisting of a massive particle 
coupled to two harmonic springs in contact with an external thermal reservoir, represented by a Langevin force term. 
External work can be done by pulling one of the springs at a given rate, which is the protocol we follow. The main 
advantage of this model is that we can explicitly make all the calculations with no approximations. No approximations 
are needed in respect to the mass of the particle, the calculations being able to take care of the particle's inertia exactly. 
This model can be thought as a controllable, and simple, heat engine. 

To the best of our knowledge, for the first a Langevin model was exactly integrated, taking into account inertia 
and general initial conditions, for a Brownian particle under the action of a given protocol. An exact form for the 
nonequilibrium work generating function (NEWGF) is obtained. The Jarzynski equality is then explicitly verified, 
such as predicted 0, showing that the method used in this paper can be seen as a first principle exact, and nontrivial, 
verification of the JE. This shows the appropriateness of using white Gaussian noise to represent the interaction 
between a thermal bath and a system. Furthermore, the work probability distribution is derived explicitly for this 
case and shows that the moments of the work W are complex functions of the work rate A and the time interval r. 
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Appendix A: Finite-time coefficients 

The exact finite time coefficients are listed below. 

k'L (2Xuj 2 -e 2 X-9) Tta+ ,jQ / rV4w 2 -6» 2 \ 
CT = . — e 2 A sin 



k'L (l + 9\) _zS2±i±i [TyJ\uj 2 -e 2 
+ (l + 9X + X 2 u 2 f 2X C0S { 2 

k'L (i + ex) 



i + ex + x 2 uj 2 



k'L (2 + eX) _x i 2 ±i A) . / rV4w 2 - 9 2 

CI = , C 2 A Sin 

2 (1 + 6X + X 2 w 2 )y/Auj 2 -9 2 \ 2 



k'L X t(h»» / t \J\uj 2 - 9 2 
e 2A cos 



aw; 



(i + ex + x 2 uj 2 ) 

k'XL„ 
1 + 9 X + X 2 oj 2 

6k' 2 L 2 (3lj 2 X 2 -6 2 X 2 + 1) T (2 +8 a) (t^/IJ 2 

e 2A sin 



(f + X 2 io 2 ) 2 - 6 2 X 2 



muj 2 y/Auj 2 — B 2 
k' 2 L 2 (-9 2 X 2 + uj 2 X 2 + 1) T( 2 +8 A) [tVAl 



(f + A 2 w 2 ) 2 - 2 A 2 



e 2 A cos 



fc' 2 i 2 (-2w 2 A 2 +uj 2 X 2 ei - 29 X- 2} _ T 
2mw 2 (l + (9A + w 2 A 2 ) ° 
fc' 2 L 2 A 2 



{sX + eX + l)k' L Q _ LS 2 ±1M t^4uj 2 -9 2 , 
^ (5) = ^ s + s 2 + . 2 )(flA + l + . 2 A 2 ) C 2A C0S ' « I (A7) 



(Al) 



(A2) 



e-— (A3) 



2m(l-(9A + u 2 A 2 ) 

k f Ln t fsA-i) 

' m (p s + + u)*) (sA — 1 j 

W ^> 2 {S) = ~ m(sX-l)(9X + l+u 2 X 2 ) ( A5) 

(-9Xs-9 2 X + 2Xuj 2 -2s-9)k' L T <2 +9 .) f tV±lu 2 -9 2 \ 

WIJs) = ) - e ^ sin (A6) 

" 3W mV4w 2 - 9 2 (9 s + s 2 + uj 2 ) {9 X + 1 + u 2 X 2 ) \ 2 J y J 
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jTL a 2 k' 2 (-9 A - 1 + Xu) (1 + 9 X + Aw) 



?ti 2 wV4w 2 -6> 2 (1 + (9 A + cj 2 A 2 ) 2 



* sin 



sin (r v^w 2 -!? 2 ) 



7 TL 2 fc /2 (-0 3 A 2 - g - 2 fl 2 A + 4 A w 2 + 3 9 X 2 lo 2 ) 
7 Ti 2 fc' 2 (6» 3 A - 4 6» A cj 2 - 4 w 2 + 4 c- T 



cos 



m 2 w 2 6» (4w 2 - 6» 2 ) (1 + 9 A + w 2 A 2 ) 



(A8) 



2 1 TL 2 k' 2 X 2 (9 2 X 2 -2lu 2 X 2 -2) 



/l2 = 



to 2 V4w 2 -6> 2 (l + (9A + w 2 A 2 ) (1 + A 2 w 2 ) 2 - 6» 2 A 2 



_ t (2 + A) 

e 2 * sin 



rV4w 2 -6» 2 



2 1 TL Q 2 k' 2 X A 9 



{l + X 2 u 2 ) 2 -9 2 X 2 
1 TL Q 2 k' 2 X* 



2 (1 + A 2 cj 2 ) 2 - 6> 2 A 2 



t (2+e a) 

e 2X cos 



e A 



rV4w 2 -6» 2 



(A9) 



^22 — 



jTL 2 k' 2 X 3 



(1 + 6»A + w 2 A 2 ) m 2 



(A10) 



fc'Uo 2 (l+w 2 A 2 -A 2 2 ) T(2+f ,A) 

■ e 2 A cos 



V4w 2 -6» 2 



(l + u; 2 A 2 ) 2 -6> 2 A 2 



k' 2 L Q 2 9(l- X 2 9 2 +3uj 2 X 2 ) 



mo; 2 v / 4a; 2 — # 2 



(l + u; 2 A 2 ) 2 -6» 2 A 2 



x (2+8 a) / r V4a; 2 — 
e ^ sin 



L 2 k'uj 2 \k'X 2 (1 + 9X) + to (-2 + (9 2 A 2 - 2A 2 w 2 - A 4 cj 4 )1 _^ 

— e A 



2muj 2 



{1 + lj 2 X 2 ) 2 -6> 2 A 2 



L 2 fc' [-A;' 2 + /c'A 2 (# 2 - A 2 w 4 ) + w 2 (to - ?ti# 2 A 2 + 2mA 2 w 2 + A 4 w 4 )] 



(l+w 2 A 2 ) 2 -6» 2 A 2 



L 2 k'uj 2 (1-9X + X 2 uj 2 ) (m + to6»A - kX 2 + mX 2 uj 2 ) 



e * 



2tow 2 (1 + cj 2 A 2 ) 2 -6I 2 A 2 



(All) 



n 2 



(1 + uo 2 X 2 - X 2 9 2 ) TL 2 k' 2 _ l2± n) 

■ e 



(1 + cj 2 A 2 ) 2 -0 2 A 2 



2 * cos 



V4w 2 -6» 2 



6Ti 2 fc' 2 (l-A 2 6i 2 + 3w 2 A 2 ) ,(2+sa) / T ^4w 2 -6» 2 ' 

■ e 2 >• sin 



mwV4w 2 - 9 2 



Tk' 2 L 2 (1 + w 2 A 2 - 9uj 2 X 3 - 2 A 2 ) 



(1 + cj 2 A 2 ) 2 -6» 2 A 2 



2tooj 2 



Tk /2 L Q 2 (1 + A 3 cj 2 + w 2 A 2 - A 2 2 ) 



(1+cj 2 A 2 ) 2 -6> 2 A 2 



2tow 2 (1 + cj 2 A 2 ) 2 -6> 2 A 2 



(A12) 
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at t < 0. In fact, the external work done by the force, described by the equation above, corresponds to the work done by 
the batteries in order to initially charge the capacitors up to a fixed voltage (part of that work shifts the particles potential 
energy and the remaining is used to create the electric field between the plates), producing a constant force upon the 
particle. It has the effect of shifting the initial potential energy of the charged particle by an amount — foxo, where xq is 
the initial position of the charged particle (if we include the batteries within the system, then the potential energy they 
provide is taken to be internal) . That external work is then given by 

rt 2 >0 

w ext,o = ~fo x = - dtfo Q(t)x t , 
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consistent with references @, 0|- The important point to be consistent with is that the external work is the one changing 
the energy landscape for the phase-space point. In contrast, if we keep the batteries as an external agent for all t, then the 
total external work is given by 

Wext = ~fa(xt 2 - xo) + Wext.o = ~f° x *%- 
The definition of what is the internal Hamiltonian, and what is external, is crucial as to which Work Fluctuation relation 
can be derived. [3, [i| 



